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We propose a quantum mechanical approach to noise in resonant tunneling 
structures, that can be applied in the whole range of transport regimes, from 
OO ■ completely coherent to completely incoherent. In both limiting cases, well 

known results which have appeared in the literature are recovered. Shot noise 
reduction due to both Pauli exclusion and Coulomb repulsion, and their com- 
' bined effect, are studied as a function of the rate of incoherent processes in the 

well (which are taken into account by means of a phenomenological relaxation 



m 
o 



time), and of temperature. Our approach allows the study of noise in a variety 
of operating conditions (i.e., equilibrium, sub-peak voltages, second resonance 
\^ \ voltages), and as a function of temperature, explaining experimental results 

^ ' and predicting interesting new results, such as the dependence of noise on 

' filled emitter states and the prediction of both increasing and decreasing shot 

^ . noise with increasing temperature, depending on the structure. It also allows 

the determination of the major contributions to shot noise suppression by 
' performing noise measurements at the second resonance voltage. 

o ■ 
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I. INTRODUCTION 

In recent years, researchers have shown increasing interest in noise in resonant tunneling 
structures. In fact, from an experimental point of view, noise measurements can provide 
information about the structure and the transport properties of resonant tunneling devices 
complementary to those given by DC characteristics and small signal AC responses. On 
the other hand, the correct prediction of noise properties is a good check for the validity of 
transport models for such devices. 

In 1989, LesovikB predicted that, in the case of completely coherent transport, shot noise 
could be lower than its classical value for totally independent electron crossings through 
the structure, i.e., the so called full shot noisep in 1990, Li and coworkers! showed the 
first experimental evidence of such phenomenon in double barrier diodes. Since then, many 
theoretical studies appeared in the literature, based on both coherentB^ and semiclassical 
models,iHil while few experimental results are availablcHnHlUli If the time of observation 
T is much longer than the average time tt an electron takes to traverse the double barrier, the 
noise spectral density at low frequencies (<C l/^r) can be reduced only if consecutive current 
pulses are correlated, i.e., if pulse distribution is sub-poissonian.i Two are the mechanisms 
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which have been considered responsible for introducing such correlation: Pauli exclusion, 
and electrostatic repulsion,0'0 which both tend to smooth fluctuations of the number of 
electrons in the well region. 

In agreement to what seems to be confirmed by experimental measurements, most the- 
oretical studies predict a maximum suppression of one half of the classical shot noise value, 
that can be obtained if the transmission probabilities of the two barriers are equal.lli How- 
ever, most astonishing appears the fact that such results have been obtained both with 
coherentBMIlZl and semiclassicali'i models, and even if time correlations between consecu- 
tive traversals of the two barriers are discarded. 

In this paper, we propose an approach addressing noise properties of generic resonant 
tunneling structures in the whole range of transport regimes, from completely coherent to 
completely incoherent, from a quantum mechanical viewpoint. In the limit of coherent 
transport the result of Lesoviklll is recovered, while, in the opposite limit of loss of coherence 
for all electrons traversing the well, the semiclassical results of Daviesi are obtained. We 
also consider the combined effects of Pauli exclusion and of Coulomb repulsion, and show 
that maximum shot noise suppression of one half is to be expected independently of the 
coherence of transport, at least up to a given amount of collisions in the well. 

We also study noise behaviour of resonant tunneling structures in various bias conditions, 
i.e., equilibrium, sub-peak voltages, and second resonance voltages. In particular, as a 
check for our model, we recover the Johnson-NyquistS noise at equilibrium. Moreover, we 
study noise dependence on temperature, which has been measured experimentally,lll but has 
received little attention from a theoretical point of view. 

The outline of the paper is as follows: in Sec. II we discuss our model for transport in 
resonant tunneling structures and for transitions through the barriers, and introduce the 
simplifications and approximations needed to address the problem analytically. In Sec. Ill 
we calculate time-dependent currents and the current power spectral density, while in Sec. IV 
we obtain noise in typical operating conditions, i.e., various applied voltages and operating 
temperatures. The particular case of large well structures (where the characteristic time for 
fluctuations of the number of carriers in the well is larger than the time of observation) is 
addressed in Sec. IV. The Conclusion section ends the paper. 



II. MODEL 

In a recent paper,@ it has been shown that the sequential tunneling approach can be 
adopted to describe the whole range of transport regimes in resonant tunneling structures, 
from purely coherent to completely incoherent. This approach is applicable when the struc- 
ture can be seen as consisting of three isolated regions fli, fl^, and fir, i-e., the left reservoir, 
the well region, and the right reservoir, respectively, that are weakly coupled through the 
two tunneling barriers 1 and 2, as is sketched in Fig. 1, so that first order perturbation 
theory is applicable. 

Let each allowed state in fli, fi^, Qr be characterized by a set of parameters ai, a^, 
ar, respectively. The density of states and the occupation factor in region fig {s = l,w,r) 
are Ps(as) and fs{ces), respectively. Following Bardeen,ii tunneling is simply treated as an 
electronic transition between levels in different regions. Given that first order perturbation 
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theory is applicable, tunneling probabilities per unit time are given by the Fermi "golden 
rule."!^ 

We also take into account the effects of elastic and inelastic collisions in the well by means 
of a single phenomeno logical relaxation time Tcow'- an electron in the well has a probability 
dt /tcoii of experiencing a collision in the infinitesimal time interval dt, and electrons emerge 
from collisions with a thermal quasi-equilibrium energy distribution and a completely ran- 
dom phase. Based on a similar model (in which a relaxation length was used instead of a 
relaxation time), a compact formula for the density of states in a quantum well has been 
obtainedS. 

It is worth noticing that in the relaxation time approximationS all collisions are effective 
in randomizing phase and relaxing energy. For simplicity, we adopt this model, and we do 
not discuss the details of energy dependence of the relaxation time. In fact, for the purpose 
of this paper, we just need to be aware of the fact that collisions affect the density of states 
by broadening and lowering the resonance peaks,0^ and affect also the occupation factor 
fw in the well, which we divide into three components,^ 

/«,(««,) = fwoictw) + fi,{oiw) + fwi(^w), (1) 

where f^o is the Fermi-Dirac occupation probability associated to the quasi-Fermi level Ej^ 
in the well, and and are the occupation factors for electrons which have come from 
the left and the right electrode, respectively, and have not undergone a collision in the well. 

Suppose that electrons are in the well: the probability that in the time interval dt 
an electron enters the well through barrier m (m = 1,2) is gm{N)dt. Following Davies et 
a/.i we call gm{N) "partial generation rate" for barrier m. The probability that in the time 
interval dt an electron escapes from the well through barrier m is rm{N)dt, where rm{N) 
is the "partial recombination rate" for barrier m. Of course, we can also define the total 
generation rate g{N) = gi{N) + g2{N), and the total relaxation rate r{N) = ri{N) + r2{N). 



A. Generation and recombination rates 

Let Hm {tti = 1,2) be the perturbation Hamiltonian due to barrier m. We start by 
considering tunneling through the first barrier: following Bardeen,E2l we obtain the matrix 
element for an electron transition from a state \ai) in Qi to a state in Q^. We indicate 
it with Mii^ = {aw\Hi\ai) . The probability per unit time z/(;^|q,„)) that an electron in Qi 
jumps into an unoccupied state |a^) in is given by the Fermi "golden rule" (to first order 
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2tt 

l^{l^\a^)) = Y I Wllw?Plfldai. (2) 



Therefore, the transition rate gi from Vti to Vt^ is obtained by integrating ^^^"^ ^ 

unoccupied states in i.e.. 



91= j V{i^\a^))Pw{'^ - fw)da^ 



-j^ I I \Miiyj\^pipyjfi{l - fu,)daidaw. (3) 
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Analogously, we can obtain the recombination rate ri: the probability per unit time 

^(\aiu)-^i) that an electron in a state {a^) in jumps in Qi is 



where we have used Mu^^ = M^^f, now ri is easily obtained by integrating z/(|a„)^i) over 
occupied states in fi^, i.e., 



We wish to point out that both gi and ri are functionals of fy^, which appears explicitly 
in and (|]), and, through the Poisson equation, affects the total Hamiltonian of the 
system, and hence both the densities of states p^i, pi, and the transition matrix elements 
Miiyj. Since we are interested in considering the effects on shot noise of Pauli exclusion and 
Coulomb repulsion associated to electrons occupying allowed states in the well, it should 
also be noticed that the effect of Pauli exclusion is accounted for in through the term 
(1 — /^), while Coulomb repulsion affects gi and ri through the term |Mi/tt,p, which depends 
on the potential energy profile modified by the charge accumulated in the well. 

The passages from (0) to @ and the above considerations can be repeated for transitions 
through barrier 2, by simply substituing pedices 1 with 2, and I with r. 



In order to study the effects on shot noise of fluctuations in the distribution of the 
occupied states in the well, we have to make some simplifying, but justified, assumptions. 

First, we assume that the occupation factors // and fr for states in the external regions 
do not fluctuate, which means that thermalization mechanisms in these regions are highly 
effective in establishing an equilibrium distribution. Then, we also assume that the effect of 
fluctations of on the potential profile is weak enough that the densities of states in the 
three regions can be considered as constant. Therefore, fluctuating terms in gi and ri are 
fw, of course, and |Mi;it,p, which depends on via the Poisson equation. 

In realistic structures, the well region contains many states, which definitely makes not 
tractable the problem of considering gi and ri as functionals of the occupation factor of each 
level in the well. Hence, we need to make a further strong assumption: that generation and 
recombination rates depend on only through the total number of electrons in the well A^, 
defined as 




(4) 




(5) 



B. Simplifying assumptions 




(6) 



In other words, we assume 



gm = 



r. 



m 



rm{N) 



for m = 1, 2. 



(7) 
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The total generation and recombination rates are 

g{N) = g^iN) + g,iN) (8) 
r{N) = n{N)+r2{N). (9) 

The eigenvalues of are positive integers, and cannot be greater than the total number 
of states in the well Nq = J p^daw Though it is not necessary nor rigorous, we will often 
consider as a number large enough to be treated as a continuous quantity, and, for 
example, will write derivatives of functions of A^. 

On one hand, it is very reasonable to state that the self-consistent energy profile, and 
therefore |Mi;^p, depend essentially on the total number of electrons in the well, and only 
to a second order on the detailed shape of the probability density distribution, i.e., on 
which particular states are actually occupied. On the other hand, we have to note that our 
assumption discards the effect of the detailed shape of the term in both and (||). We 
shall consider possible drawbacks of this last approximation later on, when discussing the 
results obtained. 

It is worth noticing that the expressions for generation and recombination rates used by 
Davies et al, i.e., ^(A^) = (A^o - N)/Te, and r{N) = N/t^ (Eq. (3.10) of Ref. g), are a 
particular case of (|) and (^). 



C. Steady state distribution of electrons in the well 

The steady state distribution po{N) of electrons in the well can be obtained by using the 
conditions that r(0) = and A^ has to be not negative along with the detailed balance on 
the ratesii 

r{N + l)po{N + 1) = g{N)po{N), (10) 
which, by induction, yields 

po(iv)=po(o)n^^^^. (11) 

m=l V / 

One can then obtain Pq{0) after imposing probability normalization, i.e.. 

No 

E Po(iV) = 1. (12) 

Af=0 



D. First order approximation of generation and recombination rates 

If the distribution of the total number of electrons in the well is narrow enough we can 
greatly improve the tractability of the problem by linearising recombination and generation 
rates. Let A^ be the number for which g{N) = r{N) and let us define AA^ = N — N. We 
develop transition rates to first order in AA^, i.e.. 
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9iN) 




AN/Tg for AA^ < giN)Tg 
for AA^ > g{N)Tg 



(13) 



r{N) 



for AA^ < -r{N)Tr 

r{N) + AN/Tr for AN > -r(iV)r^ 



(14) 



where the characteristic times r„ and for generation and recombination are defined as 
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To 



_ ^ 

dN 
dr 
dN 



N=N 



(15) 
(16) 



N=N 



In Fig. 2 a quahtative picture of the approximation made is shown. Now, by substituting 
(|1|) and (|14D into (0), we obtain 



Po(iV + l) 



g{N) - AN/Tg 



TrNo 



PoiN) r{N) + {AN + l)/Tr 



Tg L + 1 



(17) 



if we define L = AN + g{N)Tr and No = g{N){Tr + Tg). Since, according to (|T3D and ([ID 
(Fig. 2), the number of electrons in the well cannot be lower than — r{N)Tr, for which 
the recombination rate is zero, nor greater than A^ + g{N)Tg, for which the generation rate 
is zero, then po{N) is non zero only between these limits, and L runs from to A^o- From 
(p!7|) we obtain for po{N) the binomial distribution 



Po(iV) 




9 r 



(18) 



where r ^ = ^ + ^. As can be seen, (|1^) reduces to Eq. (3.13) of Ref. |, if one assumes 
the transition rates given by Eq. (3.10) of the same reference, so that, as a consequence, Nq 
reduces to A^'o, and L to A^. 

We now make the reasonable hypothesis that all the stochastic processes we are consid- 
ering are ergodic, and indicate with (a) both the expectation value and the time average of 
any quantity a. Eq. (Ill) allows us to calculate the average value of A^ and the variance: 



{N) = Y.Npo{N)=N 
var(A^) = Y,N'^P(){N) - {Nf 



Nn 



{9)r. 



(19) 
(20) 



It is worth noticing that, according to (^), the standard deviation [var(A^)]^/^ is much 
lower than the allowed range A"o of variation for A^, therefore we can reasonably assume the 
linearization of transition rates to be applicable. Moreover, given that (A^) = A^, we have 



{r) = {g) = g{N)=gm)- 



(21) 



It is useful to linearize partial recombination rates through each barrier, analogously 
to what we have done in (|13|) and (0) for g{N) and r{N). We have to define 1/Trm = 
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{dryn/dN)\^^fj and l/xg^ = {-dgm/dN)\^^^, m = 1,2. Therefore we have (gm) = gm{N) 

and {vm) = for m = 1,2. The steady-state current (i) at the device electrodes is 
given by the net transition rate through either barriers 

{i) = q{9i-ri) =q{r2-g2)- (22) 
We also define quantities that will be used in the following paragraphs: 

rfi = t;,^ + r-/, r^' = r^,' + r^f. (23) 

Of course we have 

'^g ~ '^g2 + '^gl ) '^r — '^rl + '^r2 5 (24) 

^"' = ^r' + ^2"' = ^;' + ^r'- (25) 



E. Autocorrelation function 

The autocorrelation function c^jyit) of N(t) is defined as 

CMNit) = (AiV(O)AiV(t)). (26) 
By taking the time-derivative we have 

= {AN{OMN)-r{N)]), (27) 

where we have used the rate equationi0 dN{t)/dt = g{N) - r{N). From (0), (|g), 
and (|26|), we can write dcN^/dt = —c^n/t, from which we finally have 

Civjv = ((7)re-l*l/^ (28) 

where we have used ( pOf ) and the fact that CArAr(O) = var(A^). Let us point out that r has 
the role of characteristic time of fluctuations in the number of electrons A^. 



III. NOISE 

A. Time-dependent current 

In this section, we are going to calculate the time dependent current and its noise spectral 
density in the case of a constant voltage applied between the electrodes. According to the 
Ramo-Shockley theoremll, and to the electrokinematics theoremil that generalises it to any 
system, to the electromagnetic field and to quantum mechanics,^ when an electron tunnels 
through one of the barriers, it generates a pulse in the current of the external circuitH; the 
time integral of the current pulse associated to the traversal of barrier m, (m = 1, 2) is 
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XmQ, where Am is equal to the ratio of the voltage drop across barrier m to the total applied 
voltage. Of course, we have Ai + A2 = 1- In terms of the quasi- Fermi levels of the three 
regions, we have Ai = {Efi - Ef^)/ {Efi - Ej^), and A2 = (-E/u. - -E/r)/ {Efi - Efr). Suppose 
that we observe the system in the interval (0, T), in which the current i{t) has the form 



i{t) = Xiq 



A2g 



j > 



+ 



tf) 



(29) 



where /J^ gives the shape of the current pulse due to a single generation via barrier 1 starting 
at time t^^. Traversals of the barrier and current pulse shapes are not identical, therefore 
we have to associate a different function /J^ to each pulse. What all functions fj have in 
common is the normalization to unity, and the fact that their Fourier transform Fjiuj) is 
flat and equal to unity for frequencies much smaller than the inverse of the traversal time 
of each barrier. At such frequencies the Fourier transform of the current is 



I[uj) = J i{t) exp{—iut)dt 



Xiq 



X]exp(-za;tj ) - X^exp 



+ 



exp(— zcut^^) — exp(— icut 



92 ^ 

i > 



(30) 



The integral and the sums run over all the pulses occuring in the interval of observation 
(0,T). We wish to point out that at low frequencies (here "low" means again much smaller 
than the inverse of the traversal time of the device), the power spectral density of current 
fluctuations S{ijj) is not influenced by the pulse shape, therefore we do not expect ^(ti;) to 
be dependent on the particular values of Ai and A2 (as shown in Appendix Bl). 



B. Noise spectral density 

The power spectral density Siuj) of the current is defined as 

2 



5(..) = ^(|/(^)P)-47r(z)^5(^), (31) 



where, again, T is the time of observation. By substituting ( PP] ) in ( PT| ) we obtain 16 terms 
of the type 

(EEe^pH^(^i-^f)])- (32) 

k j 

where a, /3 = gi, 92,^1, r2. We can analyze these terms following Ref. ^: if a = P the diagonal 
terms j = k are equal to unity, and their sum gives the number of a-hops occurring from 
time to T, i.e., on average, {P)T. 
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For a ^ j3, oi a = (3 and j 7^ fc, we can define hapit) as the probability per unit time 
that a a-hop occurs at time t given that a /?-hop occurred at time 0. Therefore we can write 
im) as 



{Y, f eM-^^it-tl)]Kp{t-tl)dt). 



(33) 



As can be seen, the integral in (^) is independent of time and is simply the Fourier transform 
Hajsioj) of ha/3{t). The sum over k in ( |5BD contains, on average, {P)T terms, therefore 
becomes 



(34) 



The detailed derivation of all the correlation functions Ha^jsiuj) is reported in Appendix A. 
Substitution of terms like (|34D into S{lj) finally yields, for ^ 1/r (see Appendix Bl), 



2g2 



2 / (gi + n) , {92 + ^2) 

T ^ \- - 



To 



r? 



(35) 



As expected, at frequencies smaller enough than the inverse of the transit time of electrons 
through the whole device, the relative sizes of the current pulses corresponding to the traver- 
sal of the two barriers are not relevant, therefore the dependence of S{uj) on Ai and A2 is 
cancelled out. 

We can arbitrarily decompose the net current into a component coming from the left 
and going towards the right Ii, given by the current q{gi) entering the well through barrier 1 
multiplied by the portion (r2) / (r) exiting through barrier 2, and a component coming from 
the right Ir, so that we have 



U = q- 



{9i){r2) 



Ir 



q- 



.{92) (n) 



(r) {r) 
Substitution of (|36D into (|35| ) yields (see Appendix C) 



(36) 



2q 



+ Ir 



2r^ 
T1T2 

riT2 



^2(^2), 



(37) 



In the next section we shall apply (^) and (|37|) to derive the expression of noise at several 
operating conditions. 



IV. NOISE SUPPRESSION IN DIFFERENT OPERATING CONDITIONS 



A. Thermal Equilibrium 
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1. Conductance at equilibrium 



We shall show that, at equilibrium, the noise spectral density S{uj) given by (pTl) yields 
the well known Johnson-Nyquist0 result. This is a good test for our model, and confirms 



the fact that 
noise formula 



lot noise and thermal-equilibrium noise are special forms of a more general 
fSs^ i.e., for the system considered in this paper, Eq. (|37|). 
At equilibrium, the occupation factors in each of the three regions are given by Fermi- 
Dirac statistics, and, of course, zero net average currents flow through each barrier, i.e., 
(^i) = {ri), and (^2) = (^2). 

It is reasonable to suppose that, if we apply an infinitesimal perturbation from the equi- 
librium condition, the electrons in the three regions still obey quasi-equilibrium distributions, 
i.e., the occupation factor for each region has the Fermi-Dirac form, and is associated to a 
quasi-Fermi level which tends to that in equilibrium. 

Let us calculate the conductance G at equilibrium. We need to apply a small voltage V 
between the right and the left electrode, as sketched in Fig. 3: we have qV = Efi — Efr- 
Let us define ei = Eji — Ef^, and €2 = Ef^ — Efr'. so that we also have qV = ei + 62. The 
total current is (z) = /; — /,.. From (E2]) we have 



G 



d(i) 



dV y^^ 

(l{r2) d{gi 



n) 



(r) 



dV 



+ 



qjri) d{r2 - 92) 
(r) 



dV 



v=o 



The first derivative can be written in the form 



d{gi-ri)) 



dei 



ei=0 



dei 
dV 



v=o 



(gi) dei 
ksQ dV 



(38) 



(39) 



v=o 



where the equality comes from (|B8|) , ks is the Boltzmann constant and is the temperature. 
By using (^) and the corresponding equation for transitions through the second barrier, we 
write (|38|) as 



v=o 



(40) 



V=0/ 



If we remember that at equilibrium {r2){gi) = {ri){g2), that €1 + 62 = qV, and the definition 
(p6|) of //, we finally obtain G as 



G 



qli_ 



(41) 



2. Thermal Noise 



In order to recover the Johnson-Nyquist noise, we just need to obtain a simple relation 
between ri and T2 at equilibrium. From (|15|), (|16D, and (|25|), for the calculation of ti we have 



to change the number of electrons in the well region by a small amount. At equilibrium, as 
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we said above, the effect of this operation is to shift the quasi-Fermi level in the well with 
respect to those in the left and right regions. From (^3|) and (|39|) we have 



1 



d{ri - gi) 



dei 



(ri) dei 



dei 

= dN 
N = N 



N=N 



keQ dN 



(42) 



N=N 



We know that ei + e2 = 0, because at equilibrium the left and the right electrode are at the 
same potential, therefore by comparing (^2]) and the corresponding equation for T2 we can 
finally write 



This is an interesting result: in fact, when we put it into 
noise suppression are cancelled out, and we simply obtain 

S{u) = 2q{Ii + Ir) = AGkBQ, 



(43) 

the factors responsible for 

(44) 



i.e., the Johnson- Nyquist noise (the second equality comes from (^ij) and the fact that Ii 
at equihbrium). 



B. High bias 

When the voltage V applied between the electrodes is large enough that most of the 
electrons injected in the device come from only one of the electrodes, we say the device is 
in condition of high bias. Without loss of generality, we can assume that positive voltage 
is applied to the right electrode, therefore high bias means ((72) = 0. The left-going current 
component Ir vanishes too, and the total current (i) = Ii. From (^) we can obtain the shot 
noise factor! as 

2^=/ (45) 

It is worth noticing that, from (|25|), we have 

2r2 1 

T1T2 2 

and the equal sign holds only if ti = T2. Eq. (|46|) implies that 7 > 1/2. We also notice that, 
as expected, recombination through barrier one, corresponding to injected electrons not 
contributing to the net current, reduces the shot noise suppression. Maximum suppression 
(7 = 1/2) is obtained only when ti = T2 and (ri) is zero, which means that V is high enough 
that f2/-states in the resonant energy range of the well are fully occupied or that electrons 
in the well occupy states below the left electrode conduction band edge. 
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We wish to remember that here ti and T2 take into account both the Pauh exclusion 
principle and space charge effects; therefore, (45) means that the combined effects of both 
phenomena cannot push the shot noise factor below 1/2. 

The validity of this conclusion depends on the validity of the approximations we have 
made throughout this paper. Let us recall the most relevant ones: we have used first 
order approximation for generation and recombination rates as a function of A^. Actually, 
this is not very limiting: in fact strong suppression requires strong correlation between 
electron transitions (incorrelated transitions give full shot noise); and the more transitions 
are correlated, the less the number of electrons in the well fluctuates around its average 
value, attributing validity to our first order approximation. 

The second strong assumption is that generation and recombination rates depend on 
the total number of electrons in the well, and not on the distribution of occupied levels. 
This simplification, as we said above, was necessary to make the problem tractable, but, 
on the other hand, prevents us from evaluating the possible effects on noise of "shape" 
fluctuations in the distribution of occupied states in the well. Inclusion of these effects is 
possible if one addresses a very idealized situation, for example a well with only two allowed 
levels, where simply two parameters determine the electron distribution, as Egues et a/.@ 
have done. In that case, a minimum shot noise factor of approximately 0.45 is obtained, 
when the characteristic time for transitions between the two levels in the well is equal to 
the characteristic time for generations from the emitter, and recombinations towards the 
collector. 

However, characteristic times for energy relaxation and phase randomization in the well 
in real devices are much lower than the characteristic times of escape from the well, as 
shown by the poor peak-to- valley ratio of experimental devices, compared to that predicted 
by completely coherent models. Therefore, it is very reasonable that at any moment, for a 
given number of electrons in the well, the electron distribution is practically the equilibrium 
one, in a sort of adiabatic approximation. In such case, our assumption that g and r depend 
essentially on N is practically exact: the more inelastic collisions are efficient in establishing 
an equilibrium distribution, the more we can be confident on = gm{N), and = rm{N), 
m = 1,2. 

On the basis of these considerations, we do not expect the effects we have neglected to 
play a significant role, and the shot noise factor to drop below one half. Among available 
experimental studies, only the one by BrownO exhibits a value of 7 smaller than one half, 
but the estimated accuracy of his results is not reported. 



1. One dimensional structures 

Most experimental resonant tunneling structures can be treated as one-dimensional de- 
vices: the problems are therefore simplified and many quantities of interest can be obtained 
analytically. 

In a one dimensional structure, a state in any region can be decomposed in its longitu- 
dinal component \E), its transverse component jki), and its spin component \a), i.e., for 
s = l,w,r, \as) = \Es) (8> jkx) ® |cr). Electron transitions through either barrier conserve 
spin, longitudinal energy, and transverse wave vector, therefore the problem of calculating 
generation and recombination rates can be solved just in the longitudinal direction, and the 
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results can be then integrated over transverse wave vectors and doubled to account for spin 
degeneracy (we discard, for simplicity, other spin effects). If we discard Coulomb blockade 
effects, we can use the one dimensional transition matrix elements derived in Appendix A 
of Ref. For barrier 1 we have 

\{Ei\H^\E^)\^ =\Mun,{EuE^)\^ 

= h^vi{Ei)v^{E^)Ti{E^)5{E^ - Ei), (47) 

where vi^Ei) and Vw{Eyj) are usually called attempt frequencies (because of their resemblance 
to the classical number of bounces on the barrier per second) for the states \Ei) and \Eyj), 
respectively, and Ti is the tunneling probability of barrier 1. 

Let us define p's{Es) (s = l,w,r) as the density of longitudinal states in Qg, and pt^^t) 
the density of transverse states. From (|^) and (^7[) , we have 

= Anh j dE j d\^T:Vi{E)v^{E)Ti{E)p[{E)p'^{E) x 

PT(kT)/z(i^^,kT)[l -/^(i^;,kT)]. (48) 

We assume, for simplicity, that no size effects are present in the left electrode, and that 
the longitudinal density of states pi satisfy the condition 2Trhh'i{E)p'i{E) = u{E — Ecu) , where 
u is the step function, and Ecu is the conduction band edge of the left electrode. We can 
write 

gi = 2u^{ER)Tf X 

J dEu{E - Em)p'u, j dk^prfiil - /«,), (49) 

where Ept is an arbitrary resonant energy in the well, and Tf is defined as 

^ / dEu^TiujE - E^bi)pl I dkTpTfijl - fw) 
' " ^n,{ER) J dEu{E - Em)p'u, I dk^prfiil -/«,)' ^ ^ 

i.e., is practically an average of Ti weighted on suitable couples of states for transitions from 
Qi to Qyj. Analogoulsy, for recombination we can write: 

n = 2u^{Er)T[ X 

j dEu{E - Em)pI j d^TPrfM - fi). (51) 

where T[ is defined as 

rjnr ^ / dEu^TiujE - Em)p'u, I dkTprfwjl - fi) . 
' - u^{ER)JdEu{E-EM)plIdkT:PTfM-fi)' ^ ^ 

It is worth noticing that if has a unique strong resonance for E = Er, then we have 
Tf = T^ = T,{Er). 

Now, we need to define a few quantities of interest, namely Niq, fi, and rji. We have 
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N, 



10 = J dEu{E - Em)p'^ J dkTfi, (53) 

I ^ / dEu{E - EM)p'yj I dkrrpTfifw . 

J dEu{E - Em)p'u, I dkT:pTU ' ^ ^ 
^ / dEu{E - EM)p'yj I dkrrprfw 

JdEp'^JdkT^pTU ■ ^ ^ 

iV/o has a simple interpretation as the number of electrons with longitudinal energies greater 
then Ecbi that would be in the well if the occupation factor in the well was equal to that in 
the left electrode; fi is the average of fi over occupied states in the well above the conduction 
band edge of Qi, while t]i is the ratio of the number of electrons in the well with longitudinal 
energies greater than the conduction band edge of to the total number of electrons in 
fli (therefore < r^; < 1). From (|53[-^), (^) and ( ^TD can be written as 

g, = 2v^{En)Tl{Ni, - f^N) (56) 
n = 2iy^{Ej,)T[r]iN{l- fi). (57) 

To treat the problem analytically, we assume that fi and rji are only weakly dependent 
on the number of electrons in the well, so that the major dependence of gi and ri on 
are the explicit one in Eqs. ( p6D and (|57D, and those through the tunneling probabilities Tf 
and Tf, which are due to electrostatic repulsion; the characteristic times for generation and 
recombination processes Tg^ and obtained as 



1 , , rflnTf 
— = -(^i) 



r,, dN 
1 , , rflnTf 



+ 2uUEn)Tffir]i, (58) 

N=N 

+ 2u^{ER)TM^-fi)- (59) 



N=N 



If space charge effects are not relevant, the first righthand terms of both (^8]) and (|59D 
vanish, so that for ti we have 

- = — + — = 2uUER)vi[Tffi + T[(l - /,)]. (60) 

The same passages can be repeated for the second barrier, once we simply substitute all 
pedices / with r and 1 with 2. Since we are considering the case of high bias, where /r = 
and rjr = 1, we have (r2) = 2i'w{Eji)T2N and 1/t2 = 2vwT2. By substituting these results, 
(^), and (|60|) into (^), we obtain 

, = 1 , . (61) 



{T^ + r^,[Tff] + Tai-fiW 



It is straightforward to see that, because of the term r]iT^(l — fi) added to the denominator 



of (pTD , 7 cannot be smaller than 1/2 (a closer look at ([6l| ) could also prove that 7 cannot 
be smaller than 1 — fi/2). 

Eq. (|6lD becomes more readable if we consider particular cases. First, let us suppose 



that all available states in the well are above the conduction band edge of the left electrode 
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so that we have rji = 1. This case corresponds to apphed voltages V smaller than the first 
peak voltage of the I-V characteristic. If, in addition, p'^ has a narrow peak for E = E^, 
we also have Tf = = Ti{E^), and T2 = T2{Ep), that, after substitution in (|61|), yields 

T '^Ti{ER)T2{ER)fi ^ _ ~ TfJ^g 
^ (T2(E^j)+ri(E«))2 2 ' ^ ^ 

where T^^g is the peak tunneling probability of the double barrier in the case of completely 
coherent transport, i.e., Tfo^g = 4TiT2/(Ti + T2)^. However, the presence of T^ohe does not 
mean that complete coherence is required for the applicability of (|62D: it is simply required 
that the longitudinal density of states in the well has a resonance narrow enough that all 
quantities involved in the calculation of transition rates are practically constant in the energy 
range of the resonant peak. 

Eq. (^) shows the way the product /«T^ohe affects the value of 7: maximum suppression 
(7 = 1/2) appears when both fi and T^^hc 1' i-^-' when there is large charge accumulation 
at the emitter (high //), and symmetric barrier transmission probabilities i^cohe ~ The 



devices characterized in Ref. |T8| were designed to meet these conditions near the current 
peak, where they actually exhibit a noise factor 7 1/2. 

If the fir states that can be transmitted are completely filled, we have fi = 1, and we 
recover the well known resultUffl 

rppk 

7 = 1- (63) 



2. Temperature dependence 

Recent experimental measurements of shot noise in double barrier diodes as a function of 
temperatureta have shown reduced suppression with increasing temperature. This fact has 
received little theoretical attention. A simple explanation of this effect is provided by our 
model: as temperature increases, inelastic collisions in the well become more frequent; in 
other words, the effective mean free path gets shorter, resulting in a lower and wider resonant 
peak in the density of states and a lower peak-to-valley ratio in the I-V characteristics. 

We now have to remove the hypothesis of narrow density of states p^, while keeping 
rji = 1. Let us suppose to be in the condition of high accumulation at the emitter, in other 



words, let fi = 1. From (|6TD we obtain 



In order to better explain the meaning of (0), let us consider the following situation: 
a constant voltage is applied between the electrodes, and the temperature is progressively 
raised. Suppose that the barrier dimensions are such that TxiE^) ^ T2(Er) (this is the 
case, for example, of the devices fabricated and studied by Ciambrone et al&3). For very low 
temperatures the hypothesis of narrow density of states is valid, because inelastic collisions 
are rare, therefore Eq. (|63|) is applicable, and a suppression factor of one half is expected. 
When the temperature increases, the peak of widens, and Tf and T2 start to differ 
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from Ti{Er) = T2 In particular, we would have Tf > Ti{Er) ^ T2 > T^i 
because collisions with phonons make electron in the well relax to lower energy states, so 
that generation occurs more easily at higher energies (because higher energy states in the well 
are depopulated), while recombination occurs at lower energies (because electrons occupy 
lower energy states). It is straightforward to see from (Q) that 7 depends only on the 
ratio TJ/Tf , and is minimum when that ratio is unity. As temperature increases, this ratio 
decreases, and 7 approaches unity. 



This interpretation is supported by the experimental results shown in Ref. [I8|, where, at 
temperatures up to 155 K, shot noise suppression is smoothly dependent on temperature, 
while it rapidly vanishes at higher temperatures. 

We would like to emphasize that, if Ti{Er) < T2{Er), it is possible to observe reduced 
shot noise with increasing temperature: as a matter of fact, since Tf > Ti{Er) and T2{Er) > 
TJ, at some temperature we could have Tf ^ , and have a value of 7 close to one half. 
Let us point out that the possibility of both positive and negative dependence of shot noise 
on temperature is one of the relevant prediction of our approach. 

We also want to stress the point that shot noise cannot be used as a probe for measuring 
the coherence of transport in resonant tunneling diodes. In fact, in the so-called "sequential 
regime" @0 all electrons entering the well are inelastically scattered, but the density of states 
in the well is still narrow enough that the shot noise factor 7 is not affected by the collision 
rates and is simply given by (|62D , as in the case of completely coherent transport. 

3. Second resonance 

Very small noise reduction is to be expected after the first current peak: electrons in 
the well relax towards the resonant peak of p'^, and therefore fall below the conduction 
band edge of Qi, leaving higher states in the well mostly empty, so that Pauli exclusion is 
poorly effective in preventing generation from barrier 1. Many states in the well now have 
longitudinal energies smaller then E^bi, so that we have rji appreciably smaller than unity. 
Let us also consider the case of high accumulation, i.e., fi = l: (|6TD becomes 



What happens now is that maximum suppression can only be obtained if T2 = Tfrji. If 
energy relaxation cannot be disregarded, most electrons in the well relax towards the first 
resonant peak of p^, and therefore fall below the conduction band edge of Qi, leaving higher 
states in the well empty, which means rji ^ 0. In this case, unless barrier transmission 
probabilities differ by many orders of magnitude, no noise suppression should be obtained. 
In other words, Pauli exclusion is no more effective in preventing generation from barrier 1. 
This is the case, for example of the diodes studied by Ciambrone et al.0 

As we have said in the introduction, noise characterization can provide information about 
transport in such devices complementary to that given by DC and AC characteristics. For 
example, if the density of states in the well p'^ has a second resonant level, shot noise 
measurements can provide useful insights into the coherence of transport. In fact, in the 
case of completely coherent transport, electrons do not relax to the lower resonant peak of 
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p^, and r]i is close to one, leading to a one half shot noise factor, if the barriers have equal 
transmission probability at the second resonant energy. Otherwise, even a small rate of 
incoherent processes is sufficient to relax electrons to lower levels, depopulating the second 
resonant level and leading to an rji much smaller than 1 and to a shot noise factor close to 
1. 

Moreover, looking at the dependence of shot noise on temperature, we could also de- 
termine whether Pauli exclusion or Coulomb repulsion is the dominant cause of correlation 
between different pulses: in fact, while the effectiveness of Pauli exclusion is strongly depen- 
dent on the collision rate in the well (i.e., on temperature) Coulomb repulsion (that is not 
accounted for in ([61|)) essentially depends only on the total charge accumulated in the well. 



C. Short time of observation 

Suppose that the time of observation T is much smaller than the characteristic times 
for generation and recombination through either barriers, i.e., r ^ T. In this situation, 
consecutive sub-pulses corresponding to a single electron traversing the device are separated 
by a time longer than T and do not appear to be correlated. This is easily the case, for 
example, if the well region has macroscopic dimensions, i.e., several electron diffusion lengths; 
in such a case, noise corresponding to two single barrier diodes in series is expected. 

Let us recall that (|A7| - |A1(]| ) are obtained for ujt <C 1, which is not applicable now. 
Rather, we are in the opposite limit of ujt ^ 1, where hapit) = {a), {a, (3 = ri,r2, gi, g2), 
and Hap^uj) = 2TT6{Lj){a) , an intuitive result that can be rigorously obtained from (|A6|) . By 
substituting such form of Ha/3 in ( |B1| ) and then in (^) we have 



^ = Xl{gr + n) + Xl{g2 + r2). (66) 

The full shot noise for the single barrier m (m = 1,2) in the case of uncorrelated pulses is 
given by 5*^(1^) = 2q{q{gm) + g{rm)), therefore (p6|) can be written as 



Siuj) = XlS,iiu) + XlS2iiu). (67) 

Let us recall that is the ratio of the voltage drop accross barrier m to the voltage drop 
across the whole device, (i.e., if we refer to Fig. 3, Am = ^m/gV). Near equilibrium, if Rm 
is the differential resistance of barrier m, and R = Ri + R2 is the total device resistance, we 
have Xm = Rm/R and Sm = ^kBQ/Rm (as can be obtained from (|B^)), which, substituted 
in (p7|), yields the Johnson- Nyquist result for S{uj), as expected. 

If the bias is high enough that tranport occurs in only one direction, i.e., for example, 
('^i) = (^'2) = 0, we have Si = S2 = 2q{i), with (^) = q{gi) = q{r2) and 

where the equal sign holds if Ai = A2 = 1/2, i.e., again, if the structure is symmetric. 
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V. SUMMARY 



We have shown that our approach has a definite advantage over those existing in the 
hterature, because it can be apphed in the whole range of transport regimes. We have seen 
that the results of Ref. which are valid in the case of completely coherent transport, 

are recovered, as well as those of Davies et a/.Jl which have been obtained on the basis of a 
semiclassical model, i.e., for a high rate of incoherent processes in the well. Our model has 
enabled us to explain how similar results could be obtained from quite different models, and 
even if correlation between consecutive sub-pulses was discarded. 

Moreover, we have included in our model the combined effects of Pauli exclusion and of 
Coulomb repulsion on the suppression of shot noise, and concluded that in practical devices 
a suppression in excess of one half is not to be expected. 

We have also studied shot noise in different operating conditions. At equilibrium we 
have recovered the Johnson-Nyquist noise, which is not a new result, but a good test for 
the validity of our model. At sub-peak voltages we have predicted smaller suppression of 
shot noise, due to the empty states at the cathode. At voltages higher than that for the 
first peak, except for particular cases implying strongly asymmetric barriers, we recover full 
shot noise, because time correlations between transitions into and from the well region are 
reduced due to electron thermalization. In particular, for second resonance biases, we have 
shown that the Pauli exclusion plays no role in reducing shot noise, if collisions are effective 
in establishing a thermal equilibrium energy distribution; in that case, the study of noise 
suppression in that case helps us to determine the role played by Coulomb repulsion. 

Dependence of shot noise suppression on temperature has been observed in experiments, 
and has been addressed theoretically for the first time in this paper. We have shown that 
our model simply explains the experimental results. 

In the future, numerical simulations of realistic resonant tunneling structures based on 
our model will be performed, in order to relax some of the approximations required to treat 
the problem analytically and to make a comparison with experimental results. In addition, 
we plan to perform experiments based on some new predictions of our model: in particular, 
double barrier diodes with barriers of equal transparency at the second resonance voltage 
peak, and strongly asymmetric diodes which could exhibit enhanced shot noise suppression 
with increasing temperature. 
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APPENDIX A: CORRELATION FUNCTIONS 

In this appendix we shall derive the correlation functions ha/sit) needed in Section III, 
following Davies et a/.l. Let us start by calculating the correlation function hr^r^it), {l,m = 
1,2). We need to define the conditional probability p{N,t\M,0) that electrons are in 
the well at time t, given that there were M electrons at time 0. From it, we can write the 
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conditional probability per unit time ri{t\M, 0) that an r^-transition occurs at time t, given 
that M electrons were in the well at time 0, in the form 

n(t|M,0) = ri{N)p{N,t\M,0), (Al) 

N=0 

where Nq is the total number of states in the well. 

We already introduced the probability po{M) of having M electrons in the well, therefore 
we can write 



Wo 



(r 



m)=T. rmiM)poiM), (A2) 



M=0 



from which we obtain pr^{M) , i.e., the probability that, when an r^-hop occurs, M electrons 
are left in the well, as 

yrn) 

The probability per unit time that an r;-hop occurs at time t, given that an r^-hop occurred 
at time 0, i.e., hrir^{t), is then given by 

No 

hnr^{t)= Prr.iMMt\M,0) (A4) 
From (ImD we can write (|A1|) as 



n(t|M,o)= J2 + )p{N,t\M,0) 



N=0 



= in) + ^e-l*l/^ (A5) 

Trl 

where the second equality comes from the fact that p{N, t\M, 0) is normalized to unity and 
that the number of electrons in the well relax exponentially with time constant r to the 
mean value (A^),i as we know from the rate equation dN/dt = g{N) — r{N). 
From (|14D, ( |A^jA5|) , we now write 

I No 

V™W = ^Epo(M+1)x 

Vm) M=0 



rl \ Trmi^rn) 



sums are easily evaluated if we notice that Y1n°=oPo{N)^N = and J2nLo Po{N){^N)'^ 
var(A^). For u <^1/t, the Fourier transform Hnrmi^) of ( |A6|) is 
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X , V 2r / var(iV) ' 



(A7) 



Analogously, we can obtain the correlation functions for all the other processes, i.e., for 
Lm = 1,2: 



r, . 2r / var(A^) ^ 



2r /^^ var(A^) ~ 



(A8) 
(A9) 
(AlO) 



From (|30D , by using 

gi,g2,ri,r2), we have 



APPENDIX B: USEFUL FORMULAS 
1. Calculation of S{io) 
I) and considering the fact that {P)Haf3 = {a)Hpa 



= Xi {{g,)[l + i/,„,,M] - 2(^7i)i/.,,,,H + (ri)[l + 
+ \l {{g2)[l + - 2{g2)Hr,,,,{u;) + (r2)[l + 



(Bl) 



Substitution of (|A7HXT0D in (0) and then in (|31|) yields 

2r 

+2Ai 



2r 

— var(Ar) 



{gi + n) (l 

A2 f-(^2 + ra) + -(^1 + ri) - var(iV)^) 

^ / 2r\ 2r , / 
(^2 + ra) (^1 - — j + :^var(Ar) 



But, if we put 



+xi , , 

in ( p2|) , by writing ((7) as {gi + ri + (72 + ?"2)/2, we obtain 



(A? + 2AiA2 + A2) 



+ 



(B2) 



(B3) 



2q^ ' ^ "V T-2 rf 

which reduces to ( |35D if we simply remember that Ai + A2 = 1. Note that, as expected, 
dependence of noise at low frequencies upon the relative sizes of current pulses due to 
traversal of the two barriers, is cancelled out. 

Now, from (^) and the fact that r = ri + r2, we can write 

'2(ri) 



q{92 + 



Ar2) 
^ 2(r2) 

An) 



I I, 



+ 1 Mr + /«; 



(B4) 



By substituting ( P4|) and (^) in (|35D , we straightforwardly get Eq. 
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2. Derivation of Eq. 45 



Each state ja^) (s = l,w,r) is characterized by its total energy and a set of other 
parameters f3s- Electron transitions between regions conserve energy, hence we can write 
|Mk^|2 = \M,{El,(3^,(3i)\^{E^ - El). From (|) and (|) we can write 



gi -ri 



27V 

T 



[fi{E\(3i) - UE',MWMdE' . 



(B5) 



Differentiation of ( pSD with respect to ei = Efi — Efyj, given that at equilibrium fi = 
yields 



d{gi - n] 



dei 



£1=0 



2tt 

T 



dfu 



w ) 



dE 



fw 



ei=0 



dPyjdPidE'^ , 
(B6) 



in which, according to the Fermi-Dirac statistics which holds at thermal equilibrium, we 
have 



dE 



fw 



£1=0 



/wo(f ~ fwo) 



where /cs, 0, and fy^o are the Boltzmann constant, the temperature, and Fermi-Dirac occu- 
pation factor, respectively. Substitution of (|B7|) in (|B^ ) finally yields 



djgi - n] 
dei 



£1=0 



gi 



m 
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FIGURES 




FIG. 1. A generic resonant tunneling structure consists of three isolated regions Qi, Qyj, $7; 
weakly coupled by tunneling barriers, here indicated with 1 and 2. Coupling between different 
regions has to be small enough to be treated with first order perturbation theory. 
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FIG. 2. A qualitative sketch of the total generation and recombination rates is shown (thin 
lines), along with the linearization described in Sec. IID. Nq is the total number of states in the 
well, N is the number of electrons for which g{N) = r{N), Nq is the maximum allowed excursion 
for TV in the linearized approximation. 
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FIG. 3. The longitudinal section of a resonant tunneling device is sketched; Eji, Ef^,, and Efj. 
are the quasi- Fermi levels of $7^, 0,^, respectively. A small voltage V is applied between the 
electrodes, and €i, €2 are the the partial potential energy drops across barriers 1 and 2, respectively 
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